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Herein I discuss and propose updated satellite composites of the total solar irradiance covering the period 1978-2008. The 
composites are compiled from measurements made with the three ACRIM experiments. Measurements from the NIM- 
BUS7/ERB and the ERBS/ERBE sateUite experiments are used to fill the gap from June 1989 to October 1991 between 
ACRIMl and ACRIM2 experiments. The cUmate impUcations of the alternative satellite composites are discussed by using 
a phenomenological climate model for reconstructing the total solar irradiance signature on climate during the last four 
centuries. 
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^ 1. Introduction 

A contiguous TSI database of satellite observations ex- 
tends from late 1978 to the present, covering 30 years, that 
is, almost three sunspot 11-year cycles. This database is 
comprised of the observations of seven independent experi- 
ments: NIMBUS7/ERB [Hoyt et ai, 1992], SMM/ACRIMl 
[Willson and Hudson, 1991], ERBS/ERBE [Lee III et ai, 
1995], UARS/ACRIM2 [Willson, 1994; Willson, 1997], 
SOHO/VIRGO [FrohUch et al, 1997; Crommelynck and De- 
Witte, 1997], ACRIMSAT/ACRIM3 [Willson, 2001]. There 
exists another TSI satellite record, SORCE/TIM [Kopp et al., 
2003], but it is not studied here because it started just on Febru- 
ary 2003, and it is still too short for our purpose. None of these 
independent data sets cover the entire period of observation, 
thus a composite of the database is necessary to obtain a consis- 
tent picture about the TSI variation. Herein, we use the records 
plotted in Figure 1 . 

Three TSI satellite composite are currently available: the 
ACRIM composite [Willson and Mordvinov, 2003], the PMOD 
composite [Frohlich and Lean, 1998; Frohlich, 2000, 2006] and 
the IRMB composite [Dewitte et ai, 2004], respectively: see 
Figure 2. Each composite is compiled by using different mod- 
els corresponding to different mathematical philosophies and 
different combinations of data. 

For example, one of the most prominent differences between 
ACRIM and PMOD composites is due to the different way 
of how the two teams use the NIMBUS7/ERB record to fill 
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Figure 1: TSI satellite records. ERBS/ERBE and SOHOA'IRGO TSI records 
are shifted by -5W/m- and -lOW/m^, respectively, from the 'native scale' for 
visual convenience, (units of watts/metei~ at 1 A.U.) 



the period 1 989 .5 3- 1 99 1 . 75 , the so-called ACRIM-gap between 
ACRIMl and ACRJM2 records. The consequence is that these 
two composites significantly differ from each other, in partic- 
ular about whether the minimum of the TSI during solar cy- 
cle 22-23 (1995/6) is approximately 0.45 W/m^ higher (ACRIM 
composite) or approximately at the same level (PMOD) as the 
TSI minimum during solar cycle 21-22 (1985/6). Figure 3 
shows the difference between ACRIM and PMOD. 

The difference among the TSI satellite composites has sig- 
nificant implications not only on solar physics where the cor- 
rectness of the theoretical models must be necessarily tested 
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Figure 2: ACRIM, IRMB and PMOD TSI satellite composites. 

against the actual observations, and not vice versa, but also 
on the more general global warming debate. Phenomenolog- 
ical analyses [Scafetta and West, 2007,2008] using TSI proxy, 
satellite composites and global surface temperature records of 
the past 400 years show that solar variation has been a dominant 
forcing for climate change during both the pre- and industrial 
era. According to these analyses, the sun will likely be a dom- 
inant contributor to climate change in the future. However, the 
solar contribution to the global warming during the last three 
decades remains severely uncertain due mostly to the difference 
between the TSI satelUte composites. 

The phenomenological solar signature on the global temper- 
ature is found to match quite well 400 years of data since 1600 
[Scafetta and West, 2007-2008], but such an almost continu- 
ous matching would be abruptly interrupted since 1975 if the 
PMOD composite is adopted. Instead, by adopting the ACRIM 
composite it is still possible to notice a significant correlation 
between temperature data and the reconstruction of the solar ef- 
fect on climate (see figure 6 in Scafetta and West [2007] in and 
figure 1 in Scafetta and West [2008]). Thus, a significant frac- 
tion of the +OAK warming observed from 1980 to 2007 can be 
ascribed to an increase of the solar activity if ACRIM compos- 
ite is adopted, but almost none of it would be linked to solar 
activity if PMOD composite is adopted. Evidently, if the so- 
lar contribution is uncertain, the anthropogenic contribution to 
the global wamning during the last three decades is uncertain as 
well. Hence, determining the correct TSI composite during the 
last three decades should be considered of crucial importance. 

Note that the cUmate models adopted by the Intergovernmen- 
tal Panel on Climate Change [IPCC, 2007] do not agree with the 
above phenomenological findings and predict but a minor so- 
lar contribution to climate change during the last century and, 
in particular, during the last 30 years. However, such climate 
models assume that the TSI forcing is the only solar forcing 
of cUmate and use as TSI record the one derived from the TSI 
proxy reconstructions proposed by Lean [Lean, 2000; Wang, 
2005] which are compatible with the PMOD TSI composite 
since 1978. This becomes problematic if PMOD TSI composite 
is found to be flawed. In any case, the small climate sensitiv- 



ity to solar changes predicted by the current chmate models is 
also believed to be due to the absence of several climate feed- 
back mechanisms that may be quite sensitive to solar changes, 
in addition to TSI changes alone. Some of these phenomena 
include, for example, the UV modulation of ozone concentra- 
tion that would effect the stratosphere water vapor feedback and 
the modulation of the cloud cover due to the variation of cos- 
mic ray flux which is linked to changes of the magnetic solar 
activity [Pap et al., 2004; Kirkby , 2007]. These climate mech- 
anisms are expected to magnify the influence of a solar change 
on cUmate. 

The original ACRIM composite [Willson and Mordvinov, 
2003] has been constructed by simply cahbrating the three 
ACRIM datasets and the NIMBUS7/ERB record on the base 
of direct comparison of the entire overlapping region between 
two contiguous satelhte records. This composite does not al- 
ter the actual observations as they have been pubUshed by the 
original experimental groups. However, if some degradation 
or glitches do exist in the data, this composite is flawed for 
at least two reasons: 1) the mathematical methodology used 
for merging the two contiguous satellite records, which uses 
just the average during the entire overlapping regions between 
two records, may easily give biased estimates; 2) if the NIM- 
BUS7/ERB record presents some glitches, or degradation did 
occur during the ACRIM-gap, the relative position of ACRIM 1 
and ACRIM2 is falsified. 

The IRMB composite [Dewitte et al, 2004] is constructed 
by first referring all datasets to space absolute radiometric ref- 
erences, and then the actual value for each day is obtained by 
averaging all available satellite observations for that day. Thus, 
IRMB composite adopts a statistical average approach among 
all available observations; evidently, because the daily average 
estimate is based on a small set of data (1 , 2 or in a few cases 
3 data per day), it is not statistically robust, and this may easily 
produce artificial sUps every time data from a specific record 
are missing or added. 

The PMOD composite [Frohlich and Lean, 1998; Frohlich, 
2000, 2004, 2006] is constructed by altering the published ex- 
perimental TSI satellite data every time the PMOD team claims 
that the published data are corrupted because of presumed sud- 
den gUtches due to changes in the orientation of the spacecraft 
and/or to switch-offs of the sensors, or because of some kind of 
instrumental degradation. Some TSI theoretical model predic- 
tions [Lee III e? a/., 1995; Chapman e? a/., 1996; FroohUch and 
Lean, 1998] have been heavily used by the PMOD team to iden- 
tify, correct and evaluate these presumed errors in the pubhshed 
TSI satellite records, and these models have been changed con- 
stantly during the last 10 years. 

PMOD composite is claimed to be consistent with some TSI 
theoretical proxy models [Wenzler et al., 2006; Krivova et al., 
2007]. However, differences between the model and the PMOD 
TSI composite can be easily recognized: for example, Wenzler 
et al. [2006] need to calibrate the model on the PMOD com- 
posite itself to improve the matching, and several details are 
not reproduced. Also it can not be excluded that an alterna- 
tive caUbration of the parameters of these TSI proxy models 
may better flt the ACRIM TSI satelUte composite. Evidently, if 
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Figure 3: Relative difference between ACRIM and PMOD TSI satellite com- 
posites. 

the above theoretical models and/or the corrections of the satel- 
lite records implemented by the PMOD team are found to be 
severely flawed, PMOD is flawed as well. In any case, an ap- 
parent agreement between some theoretical TSI model which 
depends on several calibration parameters and a TSI satellite 
composite does not necessarily indicate the correctness of the 
latter because in science theoretical models should be tested and 
evaluated against the actual observations, and not vice versa. 

Herein, we construct alternative TSI satellite composites us- 
ing an approach similar to that adopted by the ACRIM team, 
that is, we do not alter the pubUshed satellite data by using 
predetermined theoretical models that may bias the composite. 
However, contrary to the original ACRIM team's approach we 
use a methodology that takes into account the evident statisti- 
cal biases that are found in the published satellite records. The 
three ACRIM records are preferred and the ACRIM gap is filled 
by using the measurements from the NIMBUS7/ERB and the 
ERBS/ERBE satellite experiments. Finally, we use these alter- 
native TSI satellite composites in conjunction with a recent TSI 
proxy reconstruction proposed by Solanki's team [Krivova et 
al., 2007] to reconstruct the signature of solar change on global 
climate using a phenomenological model [Scafetta and West, 
2007, 2008]. 

2. SIMIVI/ACREVII vs. 1VI1VIBUS7/ERB 

The first step is to compose the SMM/ACRIMl and the NIM- 
BUS7/ERB records. Note that the relative accuracy, precision 
and traceability of these two databases are radically difl'erent. In 
particular, the average error of NIMBUS7/ERB measurements 
is ±0.16W/m^ while the average error of SMM/ACRIMl mea- 
surements is ±0.04 W/m^: thus SMM/ACRIMl is significantly 
more precise than NIMBUS7/ERB. Moreover, NIMBUS7/ERB 
was not able to continuously calibrate its sensor degradations 
as ACRIMl was. NIMBUS7/ERB radiometer was caUbrated 
electrically every 12 days. For the above reasons ACRIMl 
measurements are supposed to be more accurate than the NIM- 
BUS7/ERB one and, when available, they are always preferred 
to the NIMBUS7/ERB ones. 



Figure 4: Relative difference between NIMBUS7/ERB and SMM/ACRIMl 
TSI records. 

It is necessary to adopt the NIMBUS7/ERB record for 
reconstructing the TSI record during three periods: be- 
fore 17/02/1980, from 04/11/1983 to 03/05/1984, and after 
14/07/1989. To accomphsh this we evaluate the position of 
NIMBUS7/ERB relative to SMM/ACRIMl : we plot tiiis in Fig- 
ure 4. The black smooth curve is a 91 -day moving average. 

The data shown in Figure 4 have an average of 4.3W/m^. 
The original ACRIM TSI composite [Willson and Mordvi- 
nov, 2003] is constructed in a way which is equivalent to use 
the above average value to merge SMM/ACRIMl and NIM- 
BUS7/ERB record. 

However, if the differences between SMM/ACRIMl and 
the NIMBUS7/ERB records were only due to random fluctu- 
ations around an average value, the error associated to the 91- 
day moving average values had to be about +0.02-0.03W/;7i^, 
as calculated from the measurement uncertainties. Because 
the standard deviation of the smooth data shown in Figure 
4 is significantly larger, +Q.\lWlm^, the difference between 
SMM/ACRIMl and NIMBUS7/ERB measurements is not just 
due to random fluctuations, but due to biases and trends in 
the data probably due to poor sensor caUbration of NIM- 
BUS7/ERB. 

Under tiie theoretical assumption tiiat the SMM/ACRIMl 
measurements are more accurate than NIMBUS7/ERB, the 
black 91 -days moving average smooth curve shown in Figure 
4 suggests that NIMBUS7/ERB measurements could gradually 
shift during a relatively short period of time, a few months, by 
an amount that on average is about Q.lWjnp-, and in a few cases 
can also be as large as 0.5W/m^. In fact, irregular large oscil- 
lations with periods ranging from 5 to 12 months are clearly 
visible in Figure 4. 

The above finding shows that the original methodology 
adopted by the ACRIM team to merge SMM/ACRIMl and 
NIMBUS7/ERB records is hkely inappropriate because it as- 
sumes that NIMBUS7/ERB data are statistically stationary, that 
is, unbiased, while this is not what is found in the data. By not 
taking into account this problem the ACRIM team's methodol- 
ogy can introduce significant artificial slips in the TSI satellite 
composite at the chosen merging day. 
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To reduce the errors due to the above irregular large oscil- 
lations of NIMBUS7/ERB measurements we used the black 
91 -days moving average smooth curve shown in Figure 4 to 
reduce NIMBUS7/ERB record to the level of SMM/ACRIMl 
during the overlapping period and, thus, use the NIM- 
BUS7/ERB corrected record to fill all days SMM/ACRIMI 
record misses. Before 17/02/1980 NIMBUS7/ERB record 
is shifted by -4.5W/m^, while after 14/07/1989 it is shifted 
by -4.12W/m2. This means that relative to the original 
ACRIM composite, our composite will be 0.2W/m^ lower be- 
fore 17/02/1980, and OASW/m^ higher after 14/07/1989. 

However, these values do depend on the moving average win- 
dow adopted. In fact, by increasing the window the two above 
levels will approach to the average level at -4.3W/m^, which is 
the value used by the ACRIM team to merge the two records. 
Thus, the above two estimates can have an error as large as 
0.2W/m^. 

On the contrary, the PMOD team significantly alters both 
SMM/ACRIMl and NIMBUS7/ERB records before 1986: see 
Figure 3. These corrections are not justified by the data them- 
selves, but by theoretical models, which can be erroneous 
and/or may have large uncertainties. About the SMM/ACRIMl 
data PMOD team assumes that the SMM/ACRIMl record from 
1984 to 1986 significantly degradated: as Figure 3 shows, the 
position of NIMBUS7/ERB relative to SMM/ACRIMl gradu- 
ally increases from 1984 to 1986. However, this pattern can 
be caused both by a degradation of SMM/ACRIMI, as the 
PMOD team interprets, and by an increase of sensitivity of 
NIMBUS7/ERB sensors due to undetermined factors. In fact, 
we observe that from 1988 to 1989.5 the position of NIM- 
BUS7/ERB relative to SMM/ACRIMl gradually decreased of 
the same amount of the gradual increase observed from 1984 to 
1986. This suggests that N1MBUS7/ERB record of the TSI can 
gradually vary by these large amounts. 

In any case, ACRIM team has never published an update 
of their SMM/ACRIMl record and publicly disagreed with 
the PMOD team on many occasions on this issue [Willson 
and Mordvinov, 2003] because they were not able to find 
any physical explanation for this presumed degradation of the 
SMM/ACRIMl record. Herein, we beheve that the ACRIM 
team opinion cannot be just ignored and dismissed given the 
fact they are the authors of the data. The correction im- 
plemented by the PMOD team on the SMM/ACRIMI record 
should be considered hypothetical and not taken as granted. 
In any case, this correction would not alter the position of the 
TSI minimum in 1985/1986 relative to the minimum in 1996, 
which herein is a more important issue; PMOD team's correc- 
tion would only lower the TSI maximum in 198 1/1982 by about 
Q.2W/m^ relative to the TSI satellite composite after 1986. 

About the NIMBUS7/ERB record before 1980, although 
NIMBUS7/ERB trends appear to be quite uncertain, the PMOD 
team's correction of them should be considered hypothetical as 
well because they are not justified by other satellite measure- 
ments. In particular, PMOD team believes that the large NIM- 
BUS7/ERB peak occurred during the first months of 1979 (see 
Figures 1 and 2) is an artifact due to changes in the orientation 
of the spacecraft that has to be corrected. However, we observe 
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Figure 5: NIMBUS7/ERB and ERBS/ERBE TSI records during the ACRIM- 

gap. Ttie smootli curves are 91 -day moving averages calculated on the TSI 
values only for those days available in both records. 

that TSI theoretical reconstruction proposed by Solanki [Wen- 
zler et al, 2006] shows that a large TSI peak occurred during 
the first months of 1979. Look carefully at their figures 14 and 
15 where the TSI proxy reconstruction is compared with the 
PMOD composite; during the first months of 1979 there is a 
discrepancy of about \ W/m^ between the two records. This is 
a sufficient evidence for considering the PMOD team's correc- 
tions of NIMBUS7/ERB suspicious. 

In any case, the exact TSI patterns before 17/02/1980 and 
during the ACRIM-gap should be considered highly uncertain 
because they have to be derived from low quality satellite mea- 
surements. 

3. The ACRIM-gap: 15/07/1989-03/10/1991 

SMM/ACRIMl and UARS/ACRIM2 records can only be 
bridged by using two low quality satellite records: NIM- 
BUS7/ERB and ERBS/ERBE. Figure 5 shows the two records. 
Note that NIMBUS7/ERB and ERBS/ERBE present opposite 
trend. From 1990 to 1991.5 N1MBUS7/ERB record shows an 
increasing trend while ERBS/ERBE record shows a decreasing 
trend [Willson and Mordvinov, 2003]. Thus, the only two avail- 
able satellite records are not compatible with each other and at 
least one of the two is corrupted. Note that ERBS/ERBE too 
was unable to cahbrate its sensor degradations and a direct com- 
parison with the ACRIM records reveals that the discrepancy 
between local ACRIM smooth trends and the ERBS/ERBE 
smoth trends may as large as ±0.2 W/m^: the amphtude of these 
non-stationary biases is smaller than that observed in the NIM- 
BUS7/ERB measurements, but they are still significant. More- 
over, the average error of ERBS/ERBE's measurements are the 
largest among all satellite observations: ±0.26W/m^. 

The PMOD team claims that NIMBUS7/ERB record must be 
severely corrected during the ACRIM-gap. The reasons of these 
corrections can be found in the literature. Lee III et al. [1995] 
compared the NIMBUS7/ERB dataset with a TSI proxy model 
based on a multi-regression analysis of March 1985 to August 
1989 ERBS/ERBE irradiance measurements. They concluded 
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Figure 6: Relative difference between NIMBUS7/ERB original record and 
PMOD TSI satellite composite. 

that after September 1989 NIMBUS7/ERB time series appeared 
to abruptly increase by +OAW/m^ after a switch-off of NIM- 
BUS7/ERB for four days. Another +OAW/m^ upward shift ap- 
peared to occur on April 1990. Thus, this model suggests a 
two step shift correction that, once combined, would require 
the NIMBUS7/ERB record to be shifted down by O.SW/m^ at 
the end of April 1990. 

Later Chapman et al. [1996] review the finding by Lee et 
al. [1995] and concluded that on 29/09/1989 there was an up- 
ward shift of +03 1 W/m^ and on 9/05/1990 there was an upward 
shift of +Q.37W/m^: the new proposed combined two step cor- 
rection shift had to be -0.68 W/m^ after 9/05/1990. Afterward, 
Frohhch and Lean [1998] suggested a different two step shift 
correction, that is, the NIMBUS7/ERB record had to be ad- 
justed by -0.26W/m2 and -0.32W/w2 near October 1, 1989, 
and May 8, 1990, respectively. According the latter model dur- 
ing the ACRIM-gap NIMBUS7/ERB had to be shifted down by 
0.58 W/m^. 

Finally, Frohhch [2004, 2006] revised significantly his pre- 
vious model correction of NIMBUS7/ERB data. He first ac- 
knowledged that the supposed slip on May 1990 was indeed 
difficult to really identify. Then, he substituted the two step 
correction model with a new model in which there was only 
one slip on 29 September 1989 followed by a upward linear 
trend. Figure 6 shows our analysis of these new corrections in 
their latest version: on 29/09/1989 there is a supposed sUp of 
+0.47W/m2, which is significantly larger than what was previ- 
ously estimated, the upward trend is +Q.lA2Wnr^ jyear. The 
total downward shift forced on NIMBUS7/ERB record from 
1989.5 to 1992.5 is about QMWIm^, which is significantly 
larger than what was previously estimated by Frohlich himself 
and the other groups. 

From the above studies it is evident that several opinions have 
been formulated to solve the ACRIM-gap, even by the same 
research team, and they quantitatively disagree with each other. 
The above conflicting solutions indicate that it is not so certain 
how NIMBUS7/ERB should be corrected, if some corrections 
are truly needed. 

Figure 7 shows our analysis of the comparison between 
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Figure 7: 91-day moving average of the relative difference between NIM- 
BUS7/ERB original record and ERBS/ERBE TSI satellite composite during 
the ACRIM-gap. The figure shows the correction of NIMBUS7/ERB record 
implemented by PMOD team shown in Figure 6(soUd line). 

NIMBUS7/ERB and ERBS/ERBE. The 91 -day moving aver- 
age curve of the relative difference between N1MBUS7/ERB 
and ERBS/ERBE decreases until August 1989 around the time 
when SMM/ACRIMl merges with NIMBUS7/ERB at the level 
6.29W/m^, as shown in the graph. Since the beginning of 
September 1989 to the begirming of 1990 the curve rises rapidly 
by about 0.65W/;7i^. From 1990 to 1991.5 the curve rises by 
about 0.40W/wl Finally, from 1991.5 to 1993 the curve de- 
creases slightly by about O.OSW/m^. 

The total shift from 1989.5 to 1993 is about O.llWIm^. 
Note that the error related to the single measurements is about 
±0.11W/m^. Thus, the observed difference between NIM- 
BUS7/ERB and ERBS/ERBE is significant and must be inter- 
preted as due to biases in the data that are due to uncorrected 
degradation problems in the sensors or something else that bi- 
ases the TSI record. 

Figure 7 shows also the correction implemented by the 
PMOD team on NIMBUS/ERB record [FroohUch 2004, 2006]. 
It is evident that the PMOD team believes that the observed 
difference is due to uncorrected problems occurring only on 
NIMBUS7/ERB's sensors. Even so, the correction of NIM- 
BUS/ERB record implemented by the PMOD team ( 0.86W/m2 
from 1989.5 to 1992.5) appears to be overestimated at least by 
about 0.09 W/m^ because the total shift observed during the pe- 
riod is no more than about 0. 77 W/m^. The difference seems 
to be due to the fact that PMOD team did not take into ac- 
count that the real comparison must be done with the level when 
SMM/ACRIMl merges with NIMBUS7/ERB around the mid- 
dle of 1989, and the level during this period, as indicated in the 
figure, is about 6.29 W/m^. Thus, if on 29/09/1989 a jump re- 
ally occurred in the NIMBUS7/ERB record, this has to be about 
0.30W/m^, as previously estimated by Chapman et al. [1996] 
and FroohUch and Lean [1998]. Finally, the PMOD team's 
correction with a linear increase from 29/09/1989 to 1992.5 is 
poorly observed in data shown in Figure 7; it appears to be just 
a linear simplification of the complex pattern observed in the 
figure. 
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The major problem with the interpretation of the theoretical 
studies [Lee III et ai, 1995; Chapman et al, 1996; Froohlich 
and Lean, 1998; Froohlich, 2004, 2006] claiming that NIM- 
BUS7/ERB is erroneous during the ACRIM-gap is that, al- 
though these authors did notice a difference between NIM- 
BUS7/ERB and ERBS/ERBE records, they have interpreted 
such a difference as only due to a corruption of the NIM- 
BUS7/ERB record despite the fact that ERBS/ERBE too was 
unable to continuously calibrate its sensor degradations and its 
data had larger uncertainties than NIMBUS7/ERB data. In- 
deed, the increase observed in Figure 7 during the ACRIM gap 
could result from increased ERBS/ERBE degradation relative 
to NIMBUS7/ERB, a relative increase in the sensitivity of the 
NIMBUS7/ERB sensor, or both [Willson, 1997]. 

It is important to stress that in 1992 the experimental team 
responsible of NIMBUS7/ERB record [Hoyt et al., 1992] cor- 
rected all biases in the data they could find and after that ever 
come up with a physical theory for the instrument that could 
cause it to become more sensitive. The NIMBUS7/ERB cal- 
ibrations before and after the September 1989 shutdown gave 
no indication of any change in the sensitivity of the radiome- 
ter. When Lee III et al. of the ERBS team claimed there was an 
increase in NIMBUS7/ERB sensitivity, the NIMBUS7 team ex- 
amined the issue and concluded there was no internal evidence 
in the NIMBUS7/ERB record to warrant the correction that the 
latter team was proposing (personal communication with Hoyt 
in Scafetta and Willson [2009]). Perhaps the increase between 
1989 and 1991 in Figure 7 is indications of ERBS losing sensi- 
tivity rather than NIMBUS7 gaining sensitivity. 

There are several physical reasons to believe that 
ERBS/ERBE could degrade more likely than NIMBUS7/ERB 
in particular during the ACRIM-gap. For example: a) The 
N1MBUS7/ERB cavity radiometer was in a relatively high 
altitude (about 900 km) while ERBS/ERBE was in a low earth 
orbit (ca. 200 km). It is possible that ERBS would degrade 
much faster than N1MBUS7/ERB due to more atmospheric 
bombardment of its sensor. b) During the ACRIM-gap 
ERBS/ERBE was experiencing for the first time the enhanced 
solar UV radiation, which occurs during solar maxima, and this 
too may have caused a much faster degradation of the cavity 
coating of ERBS than of NIMBUS7/ERB because NIMBUS7 
already experienced such degradation during the previous solar 
maximum; c) From the spring 1990 to May/June 1991, when 
according to Figure 7 the difference between NIMBUS7/ERB 
and ERBS/ERBE increased by about 0.40W/m^, there was 
a rapid increase of cosmic ray flux, as Figure 8 shows. 
Also the latter phenomenon might have more likely affected 
ERBS/ERBE's sensors than NlMBUS7/ERB's ones, which 
already experienced a solar maximum 10 years earher. 

Moreover, the cosmic ray count is negative-correlated to TSI 
and magnetic flux, thus its minima correspond to solar activ- 
ity maxima. Figure 8 shows that the minimum around 1991.5 
was lower than the minimum around 1989.8-1990.5. This im- 
plies that according to this record the solar activity was likely 
higher around 1991.5 flian around 1989.8-1990.5. This con- 
tradicts the pattern observed in ERBS/ERBE while confirming 
NIMBUS7/ERB pattern, as Figure 5 shows. However, other 
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Figure 8: Climax cosmic ray of the University of New Hampshire (Version 4.8 
21 December 2006). Data from http://ulysses.sr.unh.edu/NeutwnMonitor/. 

solar indexes, such as the sunspot number index, present the 
opposite scenario. Thus, it is unlikely that solar proxies indexes 
can be used to solve this issue definitely. 

Thus, unless the experimental teams find a physical theory 
for explaining the divergencies observed in their own instru- 
mental measurements and solve definitely the problem, there 
exists only a statistical way to address the ACRIM-gap problem 
by using the pubhshed data themselves. This requires just the 
acknowledgment of the existence of an unresolved uncertainty 
in the TSI satellite data. This can be done by: 

1) Assuming fliat NIMBUS7/ERB is correct and 
ERBS/ERBE is erroneous; this would imply that during 
the ACRIM-gap ERBS/ERBE record degraded and should be 
shifted upward by 0.72-0.77 Ww^. 

2) Assuming that ERBS/ERBE is correct and NIM- 
BUS7/ERB is erroneous; this would imply that during the the 
ACRIM-gap NIMBUS7/ERB increased its sensitivity to TSI 
and should be shifted downward by 0.72-0.77W/ot^. 

3) Assuming that both ERBS/ERBE and NIMBUS7/ERB 
records need some corrections. 

Note that there is no objective way to implement method n. 
3 and infinitely different solutions may be proposed. For exam- 
ple, the one proposed by the PMOD team is just one proposal 
among many others. Herein we propose that all configurations 
between case n. 1 and case n. 2 may by possible, and for case 
n. 3 we just propose an average between the methods 1 and 2 
stressing that this arithmetic average should not be interpreted 
as a better physical solution to the ACRIM-gap problem. 

Figure 9 shows the three reconstructions of N1MBUS7/ERB 
record in agreement with the above three scenarios: [A] NIM- 
BUS7/ERB data are unaltered; [C] the NIMBUS7/ERB data 
are altered in such a way that their 91 -day moving average 
curve in Figure 9 matches exactly the 91 -day moving average 
curve of ERBS/ERBE shown in Figure 5; finally, in [B] the 
NIMBUS7/ERB data are altered in such a way that their 91- 
day moving average curve matches exactly the average between 
the two 91 -day moving average curves of NIMBUS7/ERB and 
ERBS/ERBE shown in Figure 5. 
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Figure 9: Reconstructions of the NIMBUS7/ERB record during the ACRIM- 

gap in agreement with three alternative scenarios. [A] data are unaltered; [C] 
the data are adapted in such a way that their smooth component matches exactly 
the smooth component of ERBS/ERBE; [B] the data are adapted in such a way 
that their smooth component matches exactly the average between the smooth 
components in [A] and [C]. 



4. SMM/ACREMl vs. UARS/ACRIM2 

To align SMM/ACRIMl and UARS/ACRIM2 records we 
proceed as follows. First, we merge NIMBUS7/ERB record 
and its two alternative records shown in Figure 9 with the 
SMM/ACRIMl record. The merging is done by uniting the 91- 
moving average mean curves at the merging day, 03/10/1991. 

Second, we use the finding shown in Figure 10. This fig- 
ure shows the overlapping period between NIMBUS7/ERB and 
UARS/ACRIM2 records. This interval is quite short and is 
made of two separated intervals during which both satellite 
measurements were interrupted for several months. Note that 
the two intervals are not ahgned: there is a difference of about 
Q.2W/np- between the two levels. Because the standard devia- 
tion of the data is about 0.26W/m^ which is significantly larger 
than the statistical error of measure 0A6W/m^, the figure indi- 
cates that the data are not statistically stationary. 

However, it is not evident which is performing poorly: 
NIMBUS7/ERB and ERBS/ERBE or UARS/ACRIM2. Be- 
cause the difference observed between NIMBUS7/ERB and 
UARS/ACRIM2 records in [A], and between the adapted NIM- 
BUS7/ERB and UARS/ACRIM2 records in [C] (where NIM- 
BUS7/ERB record is adapted to reproduce the smooth of the 
ERBS/ERBE record) are almost equal, the first impression is 
that UARS/ACRIM2 sensors experienced a downward shp be- 
tween the two intervals by about 0.2W/np-. However, because 
both NIMBUS7/ERB and ERBS/ERBE were less able to cal- 
ibrate their sensor degradation, it is still uncertain whether it 
is UARS/ACRIM2 record that has to be corrected and, if so, 
how large this corrections should be. Indeed, given the short 
time period and that both NIMBUS7/ERB and ERBS/ERBE are 
characterized by non stationary biases as large as ±0.2W/np-, 
it is possible that during 1992 the two latter records experi- 
enced a similar upward bias. Thus, here we decided to keep 
UARS/ACRIM2 record unaltered and merge the two sequences 
using the average of the relative differences during the entire 
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Figure 10: Relative difference between NIMBUS7/ERB and UARS/ACRIM2 

TSI satellite records according to the three hypotheses discussed in the ar- 
ticle. Note that NIMBUS7/ERB record has been already merged with 
SMM/ACRIMl record. Thus, the value reported in the figure refers to the posi- 
tion SMM/ACRIMl relative to UARS/ACRIM2 in the three alternative cases. 
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Figure 11: Relative difference between ACRIMSAT/ACRIM3 and 

UARS/ACRIM2 TSI satellite records (Top) and between SOHO/VIRGO and 
ACRIM2 TSI satellite records (Bottom). The latter is shifted down by IW/nf 
for visual clarity. 



overlapping period in all three cases, as shown in the Figure 
10. The error associated with this merging is about ±0.1 W/m^. 
However, if UARS/ACRIM2 record does need to be corrected, 
the global imphcation of this correction would be that the TSI 
satelUte composite before 1992.5 should be shifted downward 
by about O.lW/rrp- in all three cases. 

5. UARS/ACRIM2 vs. ACRIMSAT/ACRIM3 

The merging between UARS/ACRIM2 and ACRIM- 
SAT/ACRIM3 is done by using the information shown 
in Figure 11 that shows the relative difference between 
ACRIMSAT/ACRIM3 and UARS/ACRIM2, and for com- 
parison, the relative difference between SOHO/VIRGO and 
UARS/ACRIM2. Note that the UARS/ACRIM2 measurements 
were interrupted from 05/06/2001 to 08/16/2001. 

The latter comparison is necessary for determining the reason 
of the discrepancy observed between ACRIMSAT/ACRIM3 
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Figure 12: The three TSI satellite conposites. The composites [B] and [C] are 
shifted by -4W/m^ and -8W/m^, respectively, for visual clarity. 



and UARS/ACRIM2 which is significantly larger than the sta- 
tistical error associated with the measurements. In fact, the 
average statistical error of UARS/ACRIM2 data is O.OlW/m^, 
while the average statistical error of ACRIMSAT/ACRIM3 data 
is 0.008W/m^. The average statistical error of the relative dif- 
ference between ACRIMSAT/ACR1M3 and UARS/ACR1M2 is 
no more than 0.018W/m^. However, the data in the figure have 
a standard deviation of about O.SW/np- which is significantly 
larger than the statistical errors. Thus, the observed difference 
between ACRIMSAT/ACRIM3 and UARS/ACRIM2 records is 
not due to random fluctuations, but to non- stationary trends in 
the data. 

Because a similar pattern appears when UARS/ACR1M2 is 
compared with both SOHOA^IRGO and ACRIMSAT/ACRIM3 
records, it is likely that UARS/ACRIM2 sensors may have been 
experiencing some problem. Perhaps the annual cycle has been 
filtered off in some way. However, these problems appear to 
have significantly modified a natural variation in the TSI data 
characterized by a time scale close to 1 year. Because the differ- 
ence between ACRIMSAT/ACRIM3 and UARS/ACRIM2 ap- 
pears to present a cyclical pattern, an accurate way to merge the 
two sequences is to evaluate the average during an entire period 
of oscillation. The period from 04/05/2000 to 05/06/2001 cov- 
ers approximately one period of oscillation, and during this pe- 
riod the average difference between ACR1MSAT/ACRIM3 and 
UARS/ACRIM2 is 1. 86 W/m^: we use this value for the merg- 
ing. As the figure shows the averages during the first and the 
second half of the cycle are 1.70W/m^ and 2.02W/m^, respec- 
tively. This suggests that our merging has an uncertainty of 
about ±0.lSW/m^. 

6. Three updated ACRIM TSI composites 

The satellites records are merged and our three TSI com- 
posites are shown in Figure 12. Table 1 suimnarizes how 
SMM/ ACRIM 1 and UARS/ACRIM2 records have to be ad- 
justed to be ahgned with ACRIMSAT/ACRIM3. 

The composite [A] shows that the 1996 minimum is about 
0.67 ± 0.1 W/ot^ higher than the 1986 minimum. The compos- 



Table 1: Position in W/m^ of SMM/ACRIMI (Al) and UARS/ACR1M2 (A2) 
relative to ACR1MSAT/ACRIM3 (A3) in the three scenarios [A], [B] and [C] as 
discussed in the test. The errors are calculated by taking into account the high- 
est uncertainty due to the statistical non-stationarity of N1MBUS7/ERB and 
UARS/ACR1M2 when tiicy merge and when UARS/ACRIM2 merges ACRIM- 
SAT/ACRIM3. Note that the error of the position of SMM/ACRIMI compared 
to UARS/ACR1M2 is ±0.1 W/m^. The error associated to the non-stationarity 
of NIMBUS7/ERB during the ACRlM-gap is described by the three scenarios 
[A], [B] and [C]. 



ite [B] shows that the 1996 minimum is about 0.28 ± Q.lW/m^ 
higher than the minimum in 1986. The composite [C] shows 
that the 1996 minimum is about 0.1 1 ± O.lW/np- lower than the 
minimum in 1986. Thus, only in the eventuality that during the 
ACRIM-gap ERBS/ERBE data are uncorrupted the two solar 
minima would almost coincide, while on average the data indi- 
cate that the TSI minimum in 1996 is 0.28 ± OAW/np- higher 
than the minimum in 1985/6. 

Note that if UARS/ACRIM2 record needs to be corrected 
during its superposition with NIMBUS7/ERB, as explained 
above, the TSI 1996 minimum relative to the TSI 1986 min- 
imum would be about O.lW/m^ higher than the above three 
estimates. Thus, if this is the case, according to the satellites 
data the difference between the two minima would be about 
0.38 ± OAW/rn^. This would further stress that the TSI satelUte 
data do indicate that TSI likely increased during solar cycles 
21-23 (1980-2002). 



7. TSI proxy secular reconstructions 

It is necessary to use reconstructions of the solar activity as 
long as possible, at least one century, for determining the effect 
of solar variations on climate. The TSI record that is possi- 
ble to obtain from direct TSI satelhte measurements covers the 
period since 1978, and this period is far too short to correctly 
estimate how the Sun may have altered climate. The reason is 
because the chmate system is characterized by a slow charac- 
teristic time response to external forcing that is estimated to be 
about 8 years (which, theoretically, can be as large as 12 years) 
[Scafetta, 2008; Schwartz, 2008]. This decadal time response 
of the cUmate requires several decade long records for a cor- 
rect evaluation of an external forcing on climate. Thus, it is 
necessary to merge the TSI satellite composites with the long 
TSI secular reconstructions, which are quite uncertain because 
they are necessarily based on proxy data, and not direct TSI 
measurements. 

Long-term TSI changes over the past 400 years since the 
17th-century Maunder minimum have been reconstructed by 
several authors, for example: Hoyt and Schatten [1997], Lean 
[2000], Wang et al. [2005] and Krivova et al. [2007]. These 
TSI proxy reconstructions are based on the sunspot number 



8 







































HS1997 






























iW 






W2005 
-L2000 


















- ^ 


K2(X)7 

















1600 1650 1700 1750 1800 1850 1900 1950 2000 

year 

Figure 13: Secular TSl proxy reconstructions by Hoyt and Schatten [1997] 
(HS1997), Lean [2000] (L2000), Wang et al. [2005] (W2005) and Krivova et 
al. [2007] (K2007). K2007 has been shifted by -SW/m^ for visual conve- 
nience. 

record, the long-term trend in geomagnetic activity, the so- 
lar modulation of cosmogenic isotopes such as '"'C and "'Be 
records, and other solar related records. These observables are 
used because they are supposed to be linked to TSI variations. 
However, it is not known exactly how the TSI can be recon- 
structed from these historical records nor whether these records 
are sufficient to faithfully reconstruct TSI changes. Thus, the 
proposed TSI secular proxy reconstructions are quite differ- 
ent from each other and show different patterns, trends and 
maxima, as depicted in Figure 13. Nevertheless, they repro- 
duce similar patterns: in particular, note the minima during 
the Maunder Minimum (1645-1715) and the Dalton Minimum 
(1790-1820), and the TSI increase during the first half of the 
20''' century. 

The TSI increase during the first half of the 20''' century is 
particularly important. In fact, because the characteristic time 
response of climate to external forcing is about 8-12 years, an 
increase of TSI during the first half of the 20''' century would 
induce a warming also during the second half of the 20'* cen- 
tury, even if the TSI remains almost constant during the second 
half of the 20"^ century [Scafetta and West, 2007]. 

The four TSI proxy reconstructions shown in Figure 13 
present different trends since 1975. The TSI reconstruction by 
Hoyt and Schatten [1997] suggests that TSI increased during 
this period, as shown in our TSI satellite composites [A] and 
[B], and in the original ACRIM TSI satellite composite. How- 
ever, the other three TSI proxy reconstructions [Lean, 2000; 
Wang et al, 2005; Ki-ivova et al, 2007] suggest that TSI did 
not change on average since 1978, as shown in our TSI satel- 
lite composite [C] and in the PMOD TSI satellite composite. 
Thus, the uncertainty that we have found in composing the TSI 
satellite records appears unresolved also by using the TSI proxy 
reconstructions because different solar proxies do suggest dif- 
ferent TSI patterns as well. 

Because the TSI satellite composites refer to the actual TSI 
measurements, we propose their merging with the TSI proxy 
reconstructions for obtaining a TSI secular record. Here, we 
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Figure 14: Merging of the secular TSI proxy reconstruction by Krivova et 
al. [2007] (black) with the three TSI satellite composites proposed in Figure 
12 (grey). The TSI reconstructions [B] and [C] are shifted by -5W/»j^ and 
— lOW/m^, respectively, for visual convenience. The merging is made by shift- 
ing the oiiginal secular TSI proxy reconstruction by Krivova et al. [2007] by 
-0.5174W/m2, -0.1295W/m2 and +0.2584W/m2 in the case [A], [B] and [C] 
respectively. 

chose the most recent TSI proxy reconstruction [Krivova et al, 
2007], which has a daily resolution, and merge it to the TSI 
satellite composites in such a way that their 1980-1990 average 
coincides. Other choices and their implications by using the 
original ACRIM and the PMOD TSI satelUte composites with 
the TSI proxy reconstructions of Lean [2000] and Wang et al 
[2005] can be found in Scafetta and West [2007]. The three TSI 
merged records herein proposed are shown in Figure 14. The 
figure shows that during the last decades the TSI has been at its 
highest values since the 17''' century. 

8. Phenomenological solar signature on climate 

The phenomenological solar signature on climate can be esti- 
mated with a phenomenological energy balance model (PEBM) 
[Scafetta and West, 2007]. PEBM assumes that the climate sys- 
tem, to the lowest-order approximation, responds to an external 
radiative forcing as a simple thermodynamical system, which is 
characterized by a given relaxation time response r with a sen- 
sitivity a. The physical meaning of it is that a small anomaly 
(with respect to the TSI average value) of the solar input, mea- 
sured by A/, forces the climate to reach a new thermodynamic 
equilibrium at the asymptotic temperature value ahl (with re- 
spect to a given temperature average value). Thus, if A/(f) is a 
small variation (with respect to a fixed average) of an external 
forcing and Arj(f) is the Earth's average temperature anomaly 
induced by A/(f), ts.Ts{t) evolves in time as: 

d^T,(t) _ aAI{t) - Ar,(0 
dt T 

A model equivalent to ([ij has been used as a basic energy bal- 
ance model [North et al, 198 1 ; Douglass and Knox, 2005], but 
herein we use TSI records as a proxy forcing. 

We implement the PEBM by imposing that the global peak- 
to-trough amplitude of the 1 1 -year solar cycle signature on the 
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surface temperature is about O-IA" from 1980 to 2002, as found 
by several authors (see IPCC [2007], page 674 for details). This 
implies that the climate sensitivity Z\ \ to the 1 1-year solar cycle 
is 



Zn = 0.11 +0.02/:/Wm" 



(2) 



as found by Douglass and Clader [2002], and Scafetta and West 
[2005]. In addition the characteristic time response to external 
forcing has been phenomenologically estimated to be t = 8 + 2 
years [Scafetta, 2008; Schwartz, 2008]. Note that Scafetta 
[2008] has also found that climate is characterized by two char- 
acteristic time constants ti = 0.40 + 0.1 and T2 =8+2, with the 
latter estimate that may be a lower limit (the upper limit being 
T2 - 12 + 3 years), but a discussion about the consequences of 
this finding is left to another study [Scafetta, 2009]. 

The value of the parameter a is not calculated theoretically 
by using the TSl as a climate forcing as usually done in the tra- 
ditional climate models. The value of a is calculated by using 
the phenomenological climate sensitivity to the 11 -year solar 
cycle found in Eq. Q by means of the following equation 



a(T)=Z„^|l+|^ 



(3) 



which solves Eq. IT]. Thus, we find that the phenomenological 
climate sensitivity to TSl changes is 



a = 0.51±0.15/:/Wm" 



(4) 



With the above value of t and a Eq. JT] can be numerically 
solved by using as input the TSl records shown in Figure 14. 
The phenomenological solar signatures (PSSs) are shown in 
Figure 15 where the three PSSs are plotted since 1600 against 
a paleoclimate Northern Hemisphere temperature reconstruc- 
tion [Moberg et ai, 2005] and since 1850 against the actual in- 
strumental Northern Hemisphere surface record [Brohan et ai, 
2006]. 

The figure shows that there is a good agreement between the 
PSSs and the temperature record. The patterns between 1600 
and 1900 are well recovered. The warming during the first half 
of the 20''' century is partially recovered. Finally, since 1978 the 
output strongly depends on the TSl behavior. If the TSl recon- 
struction [A] is adopted, a significant portion of the warming, 
about 66% observed since 70s has been induced by solar varia- 
tions, while if the TSl reconstruction [C] is adopted, almost all 
warming, about 85% observed since 70s has been induced by 
factors alternative to solar variations. If the average TSl recon- 
struction [B] is adopted, at most 50% of the warming observed 
since 70s has been induced by solar variations. 

9. Conclusion 

We have reconstructed new TSl satellite composites by us- 
ing the three ACRIM records. We have shown that different 
composites are possible depending on how the ACRlM-gap 
from 1989.5 to 1992 is solved. Our three TSl composites in- 
dicate that the TSl minimum in 1996 is at least approximately 
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Figure 15: The three phenomenological solar signatures on climate (black) ob- 
tained with Eq.^forced with the three TSl records shown in Figure 14 against a 
paleoclimate Northern Hemisphere temperature reconstruction [Moberg et at, 
2005] from 1600 to 1850, and since 1850 against the actual instrumental North- 
ern Hemisphere surface record [Brohan et at, 2006] (grey). The figure shows 
temperature anomalies relative to the 1895-1905 average. 



0.30 + 0.40W/m^ higher than the TSl minimum in 1986. And 
that the two minima would approximately be located at the 
same level only in the eventuality that the TSl ERBS/ERBE 
satellite record is uncorrupted during the ACRlM-gap, a fact 
that may be not likely. 

None of the TSl satellite composites proposed by the 
ACRIM, IRMB and PMOD teams can be considered rigorously 
correct. All three teams have just adopted alternative method- 
ologies that yield to different TSl composites, but these teams 
have ignored the unresolved uncertainty in the data that yields 
to an unresolved uncertainty in the TSl composites as well. 

Note that comparison with theoretical TSl proxy models, for 
example Wenzler et ai, [2006] and Krivova et al. [2007], can- 
not be used to resolve the issue, as the PMOD team assumes, 
because: 1) In science theoretical models have to be tested 
against the observations, not vice versa; 2) The TSl proxy mod- 
els adopt a reductionistic scientific approach, that is, they as- 
sume that some given solar observable that refers to a partic- 
ular solar measure (for example measurements from magne- 
tograms or measurements of the intensity of a given frequency 
of the spectrum) can be used to faithfully reconstruct a global 
solar measure such as the TSl; 3) The TSl proxy models do 
depend on parameters that opportunely calibrated give different 
outcomes that can, eventually, fit alternative satellite compos- 
ites. 

Thus, because it is not possible to reconstruct with certainty 
the TSl behavior during the ACRlM-gap, the TSl decadal trend 
during the last three decades is unfortunately uncertain, and any 
discussion that needs to use the TSl record has to take into ac- 
count this unresolved uncertainty. 

However, because the uncertainty in the data indicate that the 
TSl minimum in 1996 is at least approximately 0.30+0.40 W/m^ 
higher than the TSl minimum in 1986, on average the satellite 
records do suggest that TSl may have increased from 1980 to 
2000. Therefore, the sun may have significantly contributed 
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to the warming observed during the last three decades, as sug- 
gested by the phenomenological energy balance model simula- 
tions herein proposed. 

Note that a recent paper by Lockwood [2008] concludes that 
even with the adoption of the original ACRIM composite, the 
sun's contribution to the global surface warming would be neg- 
ligible during the last three decades, in contrast with the find- 
ings of Scafetta and West [2007,2008] and those presented here. 
However, Lockwood' s findings derive from his evaluation of 
the characteristic time response of the chmate to solar vari- 
ation: T = 0.8 years. This value strongly differs from the 
value herein adopted of t = 8 years and recently measured 
by Scafetta [2008] and Schwartz [2008]. The problem with 
Lockwood's short time constant is that according to the climate 
physics implemented in most climate models, the characteris- 
tic time response of the chmate varies from a few months to 
several years and even decades, as Lockwood himself acknowl- 
edges in his paper (see references there). For example, the lin- 
ear upwelling/diffusion energy balance model used by Crowley 
[2000] is characterized by a time response of about r = 10 
years. In addition, Scafetta [2008] and Schwartz [2008] have 
found that climate is indeed characterized by two characteristic 
time constants, one short with a time scale of several months 
and one long with a decadal time scale. The climate processes 
with a fast response are usually responsible for the fast fluc- 
tuations seen in the data. Instead, the chmate processes with 
a slow response are those that drive the decadal and secular 
trends observed in the global temperature. This slow chmate 
response derives from the fact that the processes that regulate 
the decadal and secular variation of climate (most of all energy 
exchange with the deep ocean and changes of the albedo due 
to the melting of the glaciers and forestation and desertification 
processes) are very slow processes, and they work as powerful 
chmate feedbacks. Thus, we believe that Lockwood's analysis 
is inappropriate because it failed to take into account the chmate 
processes with a slow time response that would be responsible 
of a strong chmate response to solar changes. However, a more 
detailed discussion about this issue, which would imply also 
an update of the PEBM presented here, is left to another study 
[Scafetta, 2009]. 
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